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JE 要 : 基于 目前 主流 气候 业务 模式 系统 (BCC_CSM 1.1, ECMWF. System 5, NCEP. CFSv 2, TCC. MRI-CGCM 3), 
从 确定 性 预报 角度 评估 多 种 模式 对 青藏 高 原 前 冬 降水 的 预测 性 能 。 结 果 表 明 :多 模式 均 能 体现 出 青藏 高 原 前 冬 
(11 一 12 月 ) 降 水 空间 分 布 型 ,但 对 降水 量 级 存在 普遍 高 估 的 现象 ,其 中 ,来 自 BCC 的 模式 能 反映 出 降水 全 区 一 致 型 
和 南北 反 回 型 主 模 态 的 时 空 演 变 特征 ,EC 模式 高 佑 第 一 模 态 的 主导 作用 。 多 模式 对 降水 历史 回报 以 正 技巧 为 主 ， 
BCC 模 式 的 预报 技巧 最 优 ,TCC 仅 对 青藏 高 原 北 部 预测 较 好 。 从 可 预报 性 来 源 的 角度 分 析 发 现 , 赤 道中 东 太 平 洋 海 
温 指 数 (Nifio 3.4) 和 印度 洋 偶 极 子 (IOD ) 正 位 相对 提升 模式 预测 技巧 具有 很 好 的 指示 性 ;BCC 模 式 能 预测 出 2018 年 
高 原 前 冬 降水 异常 趋势 ,在 于 其 对 影响 青藏 高 原 降水 异常 的 关键 环流 型 具有 一 定 预测 能 
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青藏 高 原 被 称 作 地 球 “ 第 三 极 ”" 和 “亚洲 水 塔 …， 
受 东亚 季风 、 南 亚 季 风 “ 和 西风 带 的 共同 影响 。 
高 原 冬季 降水 占 年 总 量 比例 虽 小 ,但 极端 性 强 且 年 


同 影响 。Gao 等 研究 指出 ,水 汽 通 量 中 的 动力 分 
量 对 青藏 高 原 冬 季 降 雪 发 挥 重 要 作用 。Yuan 等 
研究 指出 印度 洋 偶 极 子 (IOD) 和 厄尔尼诺 -南方 涛 


际 和 季 内 变 率 大 ,在 高 原 东 划 及 巴 并 喀 拉 山 南 缘 易 
出 现 严重 雪灾 的 ,如 2018 年 冬季 ,青藏 高 原 降雪 量 级 
大 且 持 续 时 间 长 形成 历史 破 纪录 事件 ,遭遇 近 60 a 
来 最 严重 雪灾 ,对 高 寒 牧 区 畜牧 业 生 产 交通 运输 
带 来 很 大 影响 。Zhou 等 所 基于 国际 大 气 -海洋 耦合 
模式 第 5 阶段 比较 计划 (CMIP5 ) 预 估 了 北半球 未 来 
降雪 变化 ,结果 表明 2021 一 2050 年 和 2070 一 2100 
青藏 高 原 日 降雪 量 和 强 降 雪 事 件 频率 将 会 增加 ， 
这 对 青藏 高 原 地 区 未 来 雪灾 防御 提出 新 的 挑战 。 

在 天 气 尺度 上 ,青藏 高 原 东部 冬 春 降雪 的 典型 
环流 形势 有 北 痊 南 模型 . 乌 山 内 型 .阶梯 槽 型 等 “7。 
董 文 杰 等 研究 发 现 , 造 成 青藏 高 原 东 部 强 降雪 事 
件 年 际 变化 的 主要 原因 是 不 稳定 的 波 槽 气流 发 展 ， 
与 欧 亚 中 高 纬 大 型 天 气 系统 的 活动 和 环流 调整 有 
关 。Cuo 等 "将 青藏 高 原 北部 冬季 降雪 变化 归 因 于 
北大 西洋 涛 动 (NAO) .北极 涛 动 (AO) 东亚 西风 和 急 
流 (WJ) 和 El Nifio-Southern Oscillation (ENSO) 的 共 
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动 (ENSO ) 对 青藏 高 原初 冬 积 雪 具有 重要 指示 意义 ， 
Jiang 45 ^ gf 2z 3€ H TOD 异常 可 引起 印度 洋 西 部 对 
流 活 动 并 通过 激发 波 列 在 青藏 高 原 中 部 产生 气旋 
性 环流 异常 ,有 利于 引导 来 自 热带 的 水 汽 输送 至 高 
原 中 部 引起 积 雪 异常 。Shen 等 吕 进 一 步 验 证 IOD 
FU ENSO 对 青藏 高 原 前 冬 降雪 异常 的 独立 和 共同 作 
用 ,IOD 和 FNSO 正 位 相 时 通过 影响 欧 亚 南部 型 遥 
相关 型 (SEA ) 波 列 , 在 高 原 上 空 出 现 异常 上 升 运 动 
且 对 流 层 中 低层 气温 下 降 ,为 高 原 前 冬 降雪 异常 提 
供 动 力 和 热力 条 件 。 关 于 青藏 高 原 冬 季 降 水 ( 雪 ) 
方面 的 研究 ,多 侧重 针对 其 极端 特征 变化 及 机 理 归 
因 人 研究 ,但 模式 评估 或 数值 模拟 方面 的 工作 相对 

近年 来 伴随 气候 数值 模式 的 快速 发 展 , 模 式 已 
成 为 短期 气候 预测 的 有 力 工 具 之 一 。 由 中 国 发 展 
的 BCC 模 式 对 雪 盖 参数 化 等 物理 过 程 进行 改进 “， 
发 展 了 第 二 代 全 球 大 气 环 流 谱 模式 BCC_AGCM 


基金 项 目 : 国家 自然 基金 项 目 (42065003) 和 青海 省 科技 厅 基 础 研究 计划 (2021-ZJ-757) 共 同 资助 


作者 简介 : 申 红 艳 (1979-) , 女 ,高 级 工程 师 ,主要 研究 方向 为 气候 与 气候 变化 . E-mail: ishyshen@163.com 


通讯 作者 : 温 婷 婷 . E-mail: 66076253@qq.com 


1027 一 1039 页 


http: //azr.xjegi.com 


202308.00696v1 


ChinaXiv 


1028 Tog 


2.0.1 ,在 极端 气温 事件 “”、 夏 季 雨 带 的 年 代 际 变 
化 “季节 大 气 环流 异常 ”等 方面 在 中 国 区 域 展 现 
出 较 好 的 预测 性 能 ,基于 BCC_AGCM 2.2 所 建立 的 
第 二 代 月 动力 延伸 期 预测 模式 业务 系统 (DERF 
2.0) 对 我 国 月 尺度 气温 和 降水 的 预测 能 力也 明显 提 
升 *。 在 此 基础 上 ,国家 气候 中 心 进一步 建立 了 包 
含 全 球 碳 循 环 和 动态 植被 在 内 的 海 - 陆 - 气 - 冰 多 圈 
层 耦 合 的 气候 系统 模式 BCC_CSM 1.1, 并 基于 全 球 
110 km 中 等 分 辩 率 的 BCC_CSM 1.1 m 研 发 了 第 二 
代 短 期 气候 预测 模式 系统 ”, 有 效 提 升 了 东亚 夏季 
风 环 流 和 降水 预报 的 能 力 汪 2 。 欧 洲 中 期 天 气 预报 
中 心 (ECMWF) 在 上 世纪 90 年 代 初 就 建立 了 基于 持 
续 性 海 温 异 常 强 迫 .高 分 辨 率 的 气候 预测 系统 ,并 
于 2004 年 发 展 集合 成 员 多 达 51 个 的 海 气 耦 合集 合 
预测 系统 。 美 国 气候 预测 中 心 (NCEP/CPC)2011 年 
推出 第 二 代 模 式 系统 CFSv2, 包含 了 新 的 云 - 气 溶 
胶 - 辐 射 . 陆 面 .海洋 和 海 冰 过 程 以 及 新 的 海陆 气 资 
料 同 化 系统 ,是 大 气 -海洋 - 陆 面 全 耦合 的 系统 ,该 
模式 对 于 ENSO 年 际 变 率 有 较 好 的 模拟 能 力 , 对 全 
球 热带 地 区 降水 和 表面 气温 的 预测 较 上 一 版 本 有 
一 定 提 升 。 日 本 气象 厅 (JMA) 的 MRI-CGCM 模式 ， 
该 模式 由 大 气 环 流 模式 AGCM 和 海洋 环流 模式 
OGCM 耦合 而 成 。 目 前 ,基于 模式 针对 北半球 中 高 
纬 地 区 人 研究 仍然 较 少 ,尤其 对 青藏 高 原 的 模拟 性 能 
如 何 尚 未 可 知 。 

青藏 高 原 气候 严寒 ,进入 11 月 地 表 气 温 已 普遍 
降 至 零度 以 下 ,以 固态 降水 (降雪 ) 为 主 , 极 易 出 现 
极端 强 降 雪 事 件 且 降雪 不 易 融 化 造成 雪灾 ,因而 ， 
前 冬 是 高 原 雪 灾 高 发 期 ,因此 ,本 文中 青藏 高 原 前 
冬 是 指 11 一 12 月 这 个 时 段 ,利用 目前 国内 外 主流 的 
季节 气候 预测 模式 的 历史 回报 试验 数据 ,探究 多 模 
式 能 否 模 拟 再 现 青藏 高 原 前 冬 降水 的 主要 时 空 变 
化 特征 ,通过 误差 指标 评估 模式 间 的 预测 性 能 优 
F ,并 从 可 预报 性 角度 分 析 预 报 技巧 来 源 ,以 期 为 
高 原 地 区 冬季 气候 预测 及 雪灾 防御 提供 科学 参考 。 
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1 资料 与 方法 


11 资料 说 明 

采用 目前 国内 外 主流 的 气候 模式 ,包括 来 目 国 
家 气候 中 心 (BCC) 实 时 发 布 的 季节 气候 预测 模式 
(CSM 1.1) ,欧洲 中 期 天 气 预报 中 心 (ECMWT) 的 海 
气 斐 合 模 式 集合 预测 系统 (System 5) ,美国 国家 环 
境 预 报 中 心 (NCEP) 发 展 的 第 二 代 全 耦合 系统 
(CFSv 2) ,以 及 东京 气候 中 心 (CTCC) 推 出 的 大 气 - 海 
洋 环 流 耦 合 模式 (MRI-CGCM 3)( 表 1)。 统 一 选取 
各 模式 历年 9 月 起 报 的 青藏 高 原 前 冬 回 报 试验 数据 
进行 检验 ,基于 模式 原始 输出 场 ,利用 各 模式 历史 
回报 时 段 的 平均 值 计 算 模 式 预 报 的 环流 异常 场 。 
模式 场 多 年 平均 值 计 算 时 段 :1983 一 2010 年 (BCC)、 
1981—2010 íF (ECMWF) .1982 一 2010 年 (NCEP) 和 
1981 一 2010 年 (TCC)。 通 过 双 线 性 插值 方法 将 模式 
资料 从 格 点 插值 到 站 点 , 同 实况 场 进行 对 比 。 验 证 
资料 包括 中 国 气 象 局 国家 信息 中 心 的 逐 月 降水 站 点 
观测 数据 ,从 中 提取 青藏 高 原 观 测序 列 连续 的 86 个 
站 点 资料 (图 1) ,由 于 青藏 高 原 西部 观测 站 网 稀 玻 ， 
文中 主要 讨论 其 中 东部 的 检验 结果 ,资料 时 段 为 
1961 一 2020 年 ;再 分 析 资 料 来 自 美国 国家 环境 预报 
中 心 和 大 气 人 研究 中 心 (National Centers for Enviro- 
mental Prediction/National Center for Atmospheric Re- 
search, NCEP/NCAR ) 的 逐 月 位 势 高 度 场 和 海 表 温度 
场 ,水 平分 辨 率 为 2.5"x2.5" ,时 间 范 围 自 1979 一 2020 
年 。 前 冬 降 水 为 11 一 12 月 的 累计 降水 量 , 观 测 资 料 
气候 态 为 1981 一 2010 年 的 平均 值 。 
1.2 研究 方法 
1.2.1 检验 评估 方法 采用 距 平 相 关系 数 (Anoma- 
ly Correlation Coefficient, ACC) 时 间 相 关系 数 (Tem- 
poral Correlation Coefficient, TCC) 和 均 方 根 误差 
(Root Mean Square Error, RMSE) 5 三 项 评估 指标 , 
从 确定 性 预报 的 角度 客观 定量 评估 模式 的 预测 
性 能 。 


表 1 气候 模式 信息 列表 


Tab.1 The information list of climate models 


数据 来 源 名 称 预报 时 效 模式 分 辨 率 /(°) 历史 回报 时 段 
国家 气候 中 心 (BCC ) CSM 1.1 11 个 月 1x1 1983—2010 4Æ 
欧洲 中 期 天 气 预报 中 心 (ECMWT) System 5 7 个 月 1x1 1981—2010 4 
美国 国家 环境 预报 中 心 (NCEP) CFSv 2 9 个 月 1x1 1982 一 2010 年 
东京 气候 中 心 (TCC) MRI-CGCM 3 7 个 月 1x1 1979 一 2010 年 


202308.00696v1 


ChinaXiv 


80?E 


图 例 

。 气象 站 
L1 研究 范围 
ppm 


80?E 


ChinaXiv 合 作 期 刊 


7 期 申 红 艳 等 :基于 多 模式 的 青藏 高 原 前 冬 降水 预测 性 能 评估 1029 


90°E 100°E 


90°E 100°E 


注 : 底 图 采用 自然 资源 部 标准 地 图 制作 , 审 图 号 为 GS(2016)2923 号 ,对 底 图 边界 无 修改 。 


图 1 青藏 高 原 气象 观测 站 点 分 布 


Fig. 1 The meteorological stations distribution on the Tibetan Plateau 


距 平 相关 系数 (ACC) ,用 来 反映 预测 场 同 实况 
场 的 空间 相似 程度 ,计算 公式 如 下 : 


X (Ar;  Ax)X(AF,, - M) 
ACC, 2 (1) 


TESTES 
时 间 相 关系 数 (TCC ) 能 够 表征 模式 在 每 个 格 点 
的 预报 能 力 ,可 得 到 预报 技巧 空间 分 布 ,计算 公式 
如 下 : 


ACC 和 TCC 的 取 值 范围 均 在 -1~1, 越 接近 于 1 
则 表示 预报 技巧 越 高 。 

为 比较 多 套 资料 的 效果 ,文中 采用 公式 (3) 计 
算 模 式 场 相对 实况 观测 场 的 均 方 根 误差 (RMSE ) ， 
计算 如 下 : 


RMSE = (3) 


1x 
DIO =y 


式 中 : x M y 为 参与 对 比 的 两 套 资 料 。 
1.2.3 气候 指数 计算 方法 欧 亚 型 (Eurasian pat- 
tern,EU) 遥 相关 型 指数 ,反映 西欧 上 空位 势 高 度 与 
西伯 利 亚 之 间 反 相关 ,而 与 中 国 东 北 和 日 本 一 市 则 
为 正 相 关 ””。EU 指 数 计算 如 下 : 


EU- - iz (55N,20E) + lz (55°N,75E) 


-iz ( 40'N, 145 E ) 


(4) 


式 中 : Z 为 500 hPa 位 势 高 度 场 。 
1.2.2 统计 方法 文中 所 用 统计 方法 包括 趋势 分 
析 、 相 关 分 析 和 经 验 正 交 函数 分 解 (EOF) 法 *。 


2 结果 与 分 析 


2.1 多 模式 的 历史 回报 能 力 对 比分 析 
青藏 高 原 前 冬 降 水 量 自 西北 向 东南 递增 (图 2a)， 
其 西北 部 大 范 于 地 区 前 冬 降水 量 不 足 3 mm ,东南 部 
为 高 值 区 可 达 10 mm 以 上 ,是 冬季 雪灾 比较 频 发 的 
地 区 。 青 藏 高 原 前 冬 降水 均 方 差 同样 也 呈现 出 西 
北 少 -东南 多 的 分 布 格局 (图 2b) ,在 降水 量 级 大 的 
区 域 对 应 均 方差 偏 大 ,相反 量 级 小 的 区 域 均 方差 则 
较 小 ,高 值 区 位 于 横断 山脉 东西 两 侧 , 表 明 上 述 地 
区 的 前 冬 降水 年 际 变 率 较 大 ,这 一 特征 从 前 冬 降水 
历年 变化 (图 略 ) 中 也 可 得 以 印证 ,尤其 自 上 世纪 70 
年 代 中 期 至 90 年 代 末 ,年 际 振荡 频繁 日 幅 度 较 大 。 
在 此 根据 BCC_CSM 1.1、EC_System 5, NCEP_ 
CFSV 2、TCC_ MRI-CGCM 3 模式 (分 别 简 称 BCC, 
EC、NCEP、TCC) 的 历史 回报 结果 ,分 析 模 式 对 青藏 
高 原 前 冬 降水 基本 特征 的 刻画 能 力 ,包括 气候 态 和 
空间 变 率 主 模 态 特征 。 图 3 为 青藏 高 原 前 冬 多 年 平 
均 降 水 量 的 历史 回报 结果 , 同 图 2 的 降水 实况 相 比 ， 
四 个 模式 均 能 预测 出 青藏 高 原 前 冬 降 水 西北 少 - 东 
南 多 的 分 布 格局 ,但 对 降水 量 级 有 所 高 估 ,尤其 在 
降水 量 级 较 大 的 东南 部 ,模式 和 实况 的 偏差 均 较 
大 ,来 自 BCC 和 TCC 的 模式 预测 结果 均 为 正 偏差 ， 
其 中 ,BCC 在 横断 山脉 以 西 .昆仑 山脉 =- 柴 达 木 盆地 
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(a) 降水 量 


图 2 青藏 高 原 前 冬 气候 态 平均 降水 量 (a) 和 均 方差 tb) 空间 分 布 


Fig.2 The distribution of climate mean precipitation (a) and mean square deviation (b) in early winter over the Tibetan Plateau 


(a) BCC. CSM 1.1 (b) BCC CSM 1.1 
me 2 N 


-z 


340 km 


0 340 km 
- cc 


(c) EC System 5 (d) EC System 5 


>z 


0 340 km 


—— 


0 340 km 


(e) NCEP CFSv2 " (f) NCEP. CFSv 2 


人 


>z 


0 340 km 0 340 km 


(h) TCC MRI-CGCM 3 


340 km 


降水 量 /mm 


降水 偏差 /mm Cr 
10 20 30 40 50 60 70 80 90 100 0 5 1015202530354045 50 


图 3 基于 多 模式 回报 的 青藏 高 原 前 冬 气 候 态 降水 量 ( 左 ) 及 其 差 值 场 ( 右 ) 
Fig. 3 The distribution of climate mean precipitation (left) and its deviation (right) based on models hindcasting in early 


winter over the Tibetan Plateau 
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一 带 的 偏差 较 大 ,TCC 在 横断 山脉 西 侧 的 偏差 较 大 ， 
EC FU NCEP 两 种 模式 在 高 原 大 部 分 区 域 高 佑 实际 
观测 值 , 但 分 别 存 在 小 范围 的 低估 区 ,这 两 个 模式 
和 实况 的 偏差 整体 较 小 。 
降水 主 模 态 的 模拟 刻画 能 力 是 反映 模式 预测 
性 能 的 重要 方面 ,这 里 对 青 Oe ACIE 
行 EOF 分 解 ( 图 4) ,前 三 模 态 累积 方差 贡献 率 为 
51.4%( 表 2) ,经 North it" SRI 可 
代表 青藏 高 原 前 冬 降水 的 主要 特征 , 986 — Ets 73 25 
贡献 率 为 32.7% ,在 空间 上 几乎 旺 全 区 一 一 致 型 分 布 ， 
仅 高 原 东 南 侧 的 横断 山脉 出 现 小 范围 反 相 和 荷载 区 ; 
第 二 模 态 为 南北 反 向 型 ,祁连山 区 一 柴 达 木 盆 地 同 


(a) EOF1 OBS 


(b) EOF2 OBS 
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青藏 高 原 其 他 区 域 呈 反 相 分 布 特征 ;第 三 模 态 表现 
为 三 极 型 分 布 ,在 青 泸 高 原 东 部 地 区 目 北 问 南 呈 

“ 正 - 负 - 正 ? 型 分 布 ,位 于 青藏 高 原 腹地 的 三 江 源 区 
(黄河 .长江 澜沧江 ) 同 其 他 区 域 呈 反 相 分 布 。 根 
据 四 种 模式 历年 9 月 起 报 的 前 冬 降 水 ,同样 通过 
EOF 分 解 提 取 主 模 态 , 发 现 模式 可 以 体现 出 前 两 个 
主 模 态 空间 分 布 型 (图 4) ,对 于 第 一 模 态 ,TCC EC 
和 BCC 三 种 模式 同 实况 的 第 一 载荷 向 量 场 空间 相 
关系 数 显 著 通 过 0.10 显 著 性 检验 ,模式 方差 贡献 率 
在 30% 左 右 , 其 中 ,BCC 最 接近 实况 (33.7% ) ,EC 最 
Ea (42.9% ) ,明显 高 佑 第 一 模 态 的 主导 作用 ; 除 EC 
外 ,其 余 模 式 和 实况 场 的 第 一 模 态 时 间 系 数 相 关 性 


(c) EOF3 OBS 
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图 4 基于 观测 和 多 模式 回报 的 青藏 高 原 前 冬 降 水 EOF 前 三 模 态 空间 型 


Fig. 4 The first three EOF modes of TP precipitation during early winter over the Tibetan Plateau based on observation 


and multi-model hindcasting 
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表 2 观测 和 模式 回报 的 前 三 模 态 方差 贡献 率 


Tab.2 The variance contribution rates of the first three modes based on the observation and model hindcasting 


名 称 第 一 模 态 /% 第 二 模 态 /% 第 三 模 态 /9% 累积 方差 贡献 率 /% 
OBS 32.7 11.8 6.9 51.4 
BCC. CSM 1.1 33.7 21.3 10.1 65.1 
EC System 5 42.0 15.9 8.4 67.2 
NCEP. CFSV 2 28.1 14.4 9.2 51.7 
TCC. MRI-CGCM 3 38.1 19.5 10.9 68.5 


通过 0.05 显著 性 检验 ,BCC 最 高 (0.409) ,显著 通过 
0.01 信 和 度 水 平 检验 ,TCC 次 之 ,EC 最 低 ( 表 3)。 各 模 
式 对 第 二 模 态 (南北 反问 型 ) 的 空间 描述 能 力 最 佳 ， 
空间 相关 系数 显著 通过 0.05 信和 度 水 平 检验 ,其 中 
BCC、EC 和 NCEP 通 过 0.01 的 显著 性 检验 ,BCC 最 
优 , 空 间 相 关系 数 高 达 0.755。 对 于 第 三 模 态 ,模式 
的 刻画 能 力 明 显 减弱 ,三 极 型 分 布 特征 在 模式 中 均 
未 得 到 体现 ,模式 在 青藏 高 原 西 部 的 载荷 向 量 同 实 
况 一 致 ,但 在 其 东部 的 模拟 能 力 较 弱 。 模 式 中 第 二 
和 第 三 模 态 的 时 间 系 数 相关 性 均 较 弱 , 表 明 各 模式 
对 这 两 类 空间 型 年 际 变化 特征 的 刻画 能 力 较 弱 。 


此 同时 ,BCC 相 比 其 他 模式 的 均 方 根 误 差 偏 高 ,万 
其 在 青藏 高 原 南 部 和 西北 部 地 区 ,这 同 图 3 中 的 结 
果 相 对 应 ,TCC 模 式 次 之 ,EC 模式 最 小 ; 除 BCC 模 式 
外 ,其 余 模 式 的 RMSE 量 级 和 空间 分 布 都 比较 类 似 。 

由 模式 历史 回报 场 同 实况 场 的 时 间 相 关系 数 
(TCC) 分 布 可知 ( 图 6) ,模式 对 青藏 高 原 前 冬 降 水 预 
报 以 正 技 巧 为 主 ,其 中 ,来 和 目 BCC 模 式 的 显著 正 相 
关 区 履 盖 范围 较 大 ,表明 该 模式 预报 正 技巧 较 高 ， 
已 有 研究 表明 BCC 模 式 主要 对 ENSO I) E THOSE UI 
区 域 具 有 较 高 的 预报 能 力 ” ,青藏 高 原 正 好 位 于 该 
遥 相 关 沿 南 支 路 径 所 影响 到 的 区 域 ,有 助 于 提升 模 


综 上 所 述 ,BCC 模 式 对 青藏 高 原 前 冬 降 水 主 模 
态 的 模拟 性 能 最 优 , 较 准确 刻画 出 全 区 一 致 型 和 南 
北 反 向 型 的 主 模 态 特征 ,有 旦 时 间 系 数 的 演变 和 实况 
也 比较 相符 ,因此 ,BCC 能 更 好 把 握 青 藏 高 原 前 冬 
降水 主 模 态 的 时 空 演变 特征 。TCC 和 EC 模式 对 前 
两 个 主 模 态 空间 型 的 刻画 均 较 好 ,但 EC 夸大 了 第 
一 模 态 的 方差 贡献 ,CFS 则 对 全 区 一 致 型 和 三 极 型 
分 布 具有 一 定 的 模拟 能 
2.2 历史 回报 技巧 评估 

RMSE 是 通过 衡量 模式 预测 值 和 观测 值 的 偏差 
来 反映 模式 性 能 的 优 劣 ,由 多 模式 历史 回报 的 均 方 
根 误差 (RMSE) 结 果 可 知 (图 5) ,在 青藏 高 原 南部 
RMSE 值 普遍 高 于 北部 ,印证 了 上 节 得 出 的 青藏 高 
原 北 部 的 预测 技巧 高 于 南部 ,青藏 高 原 南 部 边缘 地 
带 的 RMSE 异 和 常 偏 大 ,对 应 降水 量 级 较 大 的 区 域 ;与 


式 在 该 区 域 的 预报 技巧 。EC 在 青藏 高 原 中 西部 为 
显著 正 相 关 区 ,说 明 对 这 些 区 域 具 有 较 好 的 预测 技 
巧 ,青藏 高 原 东部 的 黄河 源 区 和 西南 部 出 现 负 技 
Jj, NCEP 和 EC 的 分 布 类 似 , 正 负 技 巧 落 区 较 一 
致 ,EC 模式 的 显著 相关 区 范围 更 大 。 来 自 TCC 模 式 
的 显著 正 相 关 区 位 于 青藏 高 原 北部 ,高 原 南 部 则 号 
负 相 关 , 横 断 山脉 附近 为 显著 负 相 关 区 ,这 同 EC 模 
式 正好 相反 ,因而 在 实际 应 用 中 对 于 该 区 域 两 家 模 
式 可 进行 互补 性 参考 。 就 时 间 相 关系 数 这 一 评估 
标准 而 言 , BCC 模式 的 预报 技巧 最 优 , EC 优 于 
NECP ,而 TCC 仅 对 高 原 北部 具有 一 和 定 的 技巧 。 
印度 洋 和 太平 洋 海 温 年 际 变 率 异常 是 气候 预 
测 的 重要 信号 ,ENSO 事件 是 季节 预测 最 重要 的 可 
预报 性 来 源 , 已 有 人 研究 表明 ,BCC 二 代 模 式 对 ENSO 
预报 的 相关 技巧 与 一 代 模 式 相 比 有 显著 提升 “”。 


表 3 观测 值 与 模式 空间 型 及 其 时 间 系数 的 相关 对 照 表 


Tab. 3 The correlation coefficients between the observation and model hindcasting on patterns and time series 


名 称 第 一 模 态 第 二 模 态 
BCC_CSM 1.1 0.359" 0.755” 
EC_System 5 0.382" 0.646 ^ 

CFS NCEP CFSV 2 0.061 0.541 
TCC TCC MRI-CGCM 3 0.5107 0.443 


注 :*、** 和 *** 分 别 表示 通过 0.10、0.05 10.01 信 度 水 平 检验 。 


第 三 模 态 第 一 时 间 系 数 第 二 时 间 系 数 第 三 时 间 系 数 
0.002 0.409" 0.303 0.131 
0.106 0.279 -0.105 0.141 

0.391" 0.305" -0.103 0.317 
0.070 0.403" 0.088 0.013 


202308.00696v1 


ChinaXiv 


(a) BCC CSM 1.1 


0 


(c) NCEP CFSv 2 


>z 


均 方 根 误差 /mm 
4 10 15 20 25 


(b) EC System 5 


340 km 


LLL 


(d) TCC MRI-CGCM 3 


ChinaXiv 合 作 期 刊 


7 期 申 红 艳 等 :基于 多 模式 的 青藏 高 原 前 冬 降水 预测 性 能 评估 1033 


>z 


35 40 45 50 55 60 


图 5 多 模式 回报 的 青藏 高 原 前 冬 降水 均 方 根 误差 (RMSE ) 分 布 


Fig. 5 RMSE skill for precipitation during early winter over the Tibetan Plateau based on multi-model hindcasting 
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注 :点 状 区 域 为 达到 95% 及 以 上 置信 水 平 检验 。 
图 6 多 模式 回报 的 青藏 高 原 前 冬 降水 时 间 相 关系 数 (TCC ) 


Fig. 6 Temporal correlation coefficient (TCC) skill for precipitation during early winter over the Tibetan Plateau based 


on multi-model hindcasting 


那么 模式 预测 性 能 同 海 温 异 常 背景 是 否 存 在 关 
EK? 选用 Nitio 3.4 海 温 指数 和 印度 详 偶 极 子 (IOD ) 
指数 分 别 代 表 太 平 洋 和 印度 洋 海 温 状况 ,由 了 BCC 
FU TCC 模式 预测 的 ACC 分 别 同 秋 季 (9 一 10 月 ) 
Niño 3.4 和 IOD 指数 的 散 点 关系 图 可 以 发 现 (图 7)， 


当 Niño 3.4 和 IOD 为 正 位 相 时 ,对 应 模式 预测 结 

以 正 技巧 为 主 ,其 中 BCC 模式 表 现 最 优 ,相反 二 者 
为 负 位 相 时 , 则 模式 出 现 正 技巧 和 负 技 巧 的 概率 相 
当 , 表 明 模 式 预 测 的 不 确定 性 增 大 ,尤其 当 IOD 为 负 
位 相 时 ,TCC 模 式 的 ACC 以 负 值 为 主 , 说 明 此 时 模 
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(b) BCC IOD & ACC 
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图 7 BCC 和 TCC 模 式 的 距 平 相关 系数 (ACC) 同 Niao 3.4 IOD 指数 的 散 点 关系 
Fig. 7 Scatter plot of the Niño 3.4, IOD indices and ACC in early winter during 1961-2018 by BCC and TCC models 


式 已 不 具备 预测 能 力 , 可见 Niño 3.4 40 IOD 的 位 相 
对 BCC、TCC 模 式 的 预报 技巧 具有 指示 作用 ,二 者 正 
位 相 均 有 利于 这 两 种 模式 呈现 出 正 技巧 ,而 负 位 相 
时 预测 不 确定 性 增 大 ,预测 性 能 随 之 降低 , H. TCC 
模式 更 依赖 于 IOD 位 相 变 化 ,IOD 负 位 相 时 该 模式 
易 出 现 负 技巧 。 
2.3 历史 典型 年 检验 评估 

青藏 高 原 2018 年 前 冬 降 水 突破 有 历史 记录 以 
来 的 历史 极 值 ,造成 严重 雪灾 ,由 2018 年 前 冬 降 水 
量 及 其 距 平 百分率 可 知 (图 8) ,降水 高 值 区 位 于 青 
藏 高 原 东南 部 , 且 大 部 分 地 区 降水 较 气 候 值 偏 多 1 
倍 以 上 ,模式 对 如 此 典型 的 降水 异常 年 的 模拟 能 


(a) 降水 量 


降水 量 /mm L—I IEEE T CNN 
5 10 15 20 25 30 35 40 45 50 


如 何 ? 对 2018 年 历史 典型 异常 年 ,检验 多 模式 对 极 
端 降水 的 预测 效果 ,进而 评估 模式 对 极端 性 的 预测 
能 力 。 

从 2018 年 10 月 模式 起 报 的 青藏 高 原 前 冬 降水 
量 及 距 平 百分率 分 布 可 知 (图 9),BCC ECHI NCEP 
模式 存在 对 青藏 高 原 西 北部 降水 一 致 高 估 东南 部 
低估 现象 ,而 TCC 模式 则 对 整个 青藏 高 原 地 区 高 
fh , 且 在 西南 侧 高 估 最 明显 。 比 较 来 看 ,BCC 同 实 
况 最 接近 ,偏差 可 控制 在 10% 以 内 ,对 降水 分 布 型 
及 大 值 中 心 具 有 预测 能 力 ;TCC 次 之 ,在 青藏 高 原 西 
南 侧 偏差 较 大 (高 于 20% ) ,其 余地 区 的 偏差 在 10% 
左右 ;EC 和 NCEP 同 实测 的 偏差 比较 明显 。 根 据 降 


(b) 降水 距 平 百分率 
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图 8 2018 年 青藏 高 原 前 冬 降 水 量 和 降水 距 平 百分率 


Fig. 8 The precipitation and its anomalous percentage during early winter 2018 over the Tibetan Plateau 
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图 9 多 模式 预测 的 2018 年 前 冬 降水 量 ( 左 ) 及 其 降水 距 平 百分率 ( 右 ) 分 布 


Fig. 9 The precipitation (left) and its anomalous percentage (right) in early winter 2018 based on multi-model prediction 


水 距 平 百分率 分 布 , BCC 模式 预 测 结果 除 在 青藏 高 
原 东 侧 和 南 侧 与 实况 相反 ,其余 大 部 分 地 区 同 实况 
场 相 一 致 ,可 预测 出 降水 异常 偏 多 ,但 其 异常 程度 
同 实况 存在 差异 ;EC 和 NCEP 仅 在 青藏 高 原 西北 诗 
预测 出 偏 多 的 特征 ,而 在 2018 年 雪灾 异常 严重 的 高 
原 东南 部 出 现 预报 失误 ;TCC 模式 结果 显示 在 青藏 
高 原 大 部 地 区 前 冬 降水 偏 多 为 主 , 但 异常 量 级 远 不 
及 实况 。 综 合 而 言 ,BCC 和 TCC 模 式 可 较 好 预测 出 
2018 年 青藏 高 原 降 水 偏 多 的 趋势 ,但 其 异常 量 级 和 
落 区 同 实况 存在 一 定 差异 ,EC 和 NCEP 模 式 对 青藏 
高 原 西 北部 偏 多 趋势 预测 准确 ,但 在 其 东南 部 出 现 
明显 误 判 。 


环流 是 直接 影响 降水 异常 的 主要 因素 ,2018 年 
前 冬 500 hPa 位 势 高 度 距 平 场 显示 (图 10a), 自 北大 
西洋 北部 -欧洲 西部 -亚洲 区 域 呈 “ 负 - 正 - 负 ” 波 列 式 
分 布 , 欧 洲 高 压 疹 明显 偏 强 , 在 斯 堪 的 纳 维 亚 上 空 出 
现 明 显 正高 度 异 常 中 心 , 同 时 在 东北 亚 地 区 存在 一 
个 正 异常 中 心 , 有 利于 鄂 瞧 次 克海 阻塞 高 压 的 发 展 ， 
而 极地 一 西伯 利 亚 为 负 异 党 并 癌 南 延伸 至 我 国 西南 
地 区 ,东亚 大 覃 偶 强 ,青藏 高 原 位 于 负 高 度 距 平 控制 
区 ,这 种 水 平 环流 异常 特征 正好 对 应 欧 亚 遥 相 关 型 
(ED) 的 负 位 相 特 征 ”。 融 相关 指数 计算 结果 显示 
(图 10b) ,2018 年 前 冬 EU 指数 为 -129.4, 为 近 58 a 
来 历史 第 三 低 值 , 仅 次 于 2011 年 (151.9) 和 1968 年 
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图 10 基于 NCEP 再 分 析 资 料 的 2018 年 前 冬 500 hPa 位 势 高 度 距 平 场 (a) 和 前 冬 欧 亚 肥 相关 型 指数 (EU ) 的 演变 曲线 (b) 
Fig. 10 The anomalous field of 500 hPa geopotential height in early winter 2018 (a) and the temporal variation of EU index based 


on NCEP dataset (b) 
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图 11 多 模式 预测 的 2018 年 前 冬 500 hPa 位 势 高 度 距 平 场 
Fig. 11 The anomalous field of 500 hPa geopotential height in early winter 2018 based on multi-model prediction 
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(-129.1) ,可 见 2018 年 前 冬 EU 型 遥 相 关 负 位 相 特征 
异常 显著 。 分 析 不 同 模式 对 2018 年 前 冬 环流 场 的 预 
il 8E 77 ,在 此 利用 NCEP 再 分 析 资 料 进行 对 比 ,BCC 
模式 可 预测 出 中 纬度 的 波 列 结构 (图 11) ,从 北大 西 
洋 至 东亚 地 区 的 “ 负 - 正 - 负 ” 波 列 式 分 布 ,大 西洋 东 
北部 为 负 异 常 斯 堪 的 纳 维 亚 半岛 一 东欧 平原 的 正 
异常 .贝加尔 湖 及 其 以 南 的 负 异 常 , 高 技巧 区 与 欧 
亚 遥 相关 型 (EU) 波 列 分 布 类 似 。BCC 模式 预测 的 
北半球 高 度 距 平 场 同 实况 场 相似 系数 为 0.533 , 显著 
通过 0.01 信 度 水 平 检 验 , 同 实况 较 吻 合 , 对 提升 降 
水 预测 性 能 具有 直接 作用 ,但 对 巴尔 喀什 湖 一 贝 加 
尔 湖 地 区 的 预测 技巧 为 负 , 对 冬季 高 原 高 度 场 异常 
会 产生 直接 影响 ,从 而 削弱 对 降水 的 预报 性 能 ,而 
其 他 模式 对 中 纬度 的 波 列 也 有 所 体现 ,但 强度 和 中 
心 位 置 存在 明显 差异 ,在 该 区 域 的 预测 性 能 也 为 负 
或 明显 减弱 ,进一步 说 明 改 善 模 式 对 冬季 中 纬度 地 
区 的 预测 性 能 至 关 重 要 ,同时 也 为 重点 改进 区 域 提 
供 一 定 思 路 。 

2018 年 海 温 异 常 监测 显示 (图 12) ,赤道 中 东 太 
平 洋 海 温 异常 偏 暧 , 白 令 海 附近 为 显著 正 异 常 , 北 
大 西洋 地 区 为 北大 西洋 三 极 子 (NAT) 正 位 相 特 征 ， 
北 印 度 洋 一 致 偏 暧 , 南 印 度 洋 具 有 西 正 东 负 型 分 布 
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特征 ;多 模式 较 一 致 地 预测 出 赤道 太平 洋 地 区 海 温 
正 异 常 的 特征 (图 13) ,同时 对 于 北大 西洋 海 温 三 极 
型 (NAT) 正 位 相 结构 也 均 有 体现 ,尤其 是 BCC 和 下 C 
模式 对 关键 区 海 温 的 模拟 同 实况 较 接近 ,以 BCC 模 
式 表现 最 突出 ,表现 在 对 南 印度 洋 偶 极 型 海 温 、 北 
大 西洋 海 温 .ENSO 暖 位 相 特征 均 和 实况 最 接近 , 综 
合 上 述 分 析 结 果 ,BCC 模 式 对 2018 年 青藏 高 原 前 冬 
降水 和 中 高 纬 波 列 结构 的 预测 技巧 均 比 较 高 ,由 此 
推测 模式 可 通过 提高 对 关键 区 海 温 的 预报 技巧 并 
准确 把 握 海 温 对 热带 以 外 地 区 的 遥 相 关 影 响 , 从 而 
提升 整体 预测 性 能 。 
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图 12 2018 年 前 冬 海 表 温 度 距 平 场 
Fig. 12 The SSTA in early winter 2018 
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图 13 2018 年 前 冬 海 温 预报 场 
Fig. 13 The SSTA prediction in early winter 2018 from models 
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3 结论 


采用 国内 外 主流 推广 的 4 种 业务 模式 ,评估 多 
模式 对 青藏 高 原 前 冬 降 水 变化 特征 及 其 异常 的 预 
测 性 能 。 主 要 结论 如 下 : 

(1) 多 模式 能 体现 出 降水 西北 多 -东南 少 的 整 
体 分 布 格局 ,但 一 致 高 估 了 降水 量 级 ;多 模式 对 全 
区 一 致 型 降水 主 模 态 的 模拟 能 力 较 好 ,来 自 BCC 的 
模式 可 较 好 模拟 出 降水 主 模 态 的 时 空 变化 特征 。 

(2) 多 模式 对 低 纬 热带 环流 的 预报 技巧 一 致 较 
好 ,而 在 中 纬度 地 区 ,各 模式 的 预报 技巧 较 低 ,以 负 
技巧 为 主 , 夏 秋季 出 现 ENSO 或 印度 洋 海 温 正 位 相 
时 模式 的 预报 技巧 会 显著 提高 , 且 预 报 技巧 对 IOD 
指数 的 依赖 性 更 高 ,初秋 IOD 正 异常 有 利于 提升 
BCC 和 TCC 模 式 的 预测 性 能 。 

(3) 对 于 2018 年 历史 典型 异常 年 ,BCC 模 式 预 
测 结果 更 接近 实况 ,偏差 可 控制 在 10% 以 内 ,是 由 
于 BCC 模 式 对 2018 年 中 高 纬 波 列 结构 的 预测 技巧 
较 高 ,对 关键 区 海 温 的 预测 同 实 况 较 接近 ,由 此 推 
测 模 式 是 通过 提高 对 关键 区 海 温 的 预报 技巧 并 准 
确 把 握 海 温 对 热带 外 地 区 的 遥 相 关 影 响 ,从 而 提升 
整体 预测 性 能 。 
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Evaluation of multi-model precipitation simulation over the 
Tibetan Plateau in early winter 


SHEN Hongyan, WEN Tingting, ZHAO Xianrong, FENG Xiaoli 
(1. Shaanxi Meteorological Bureau, Key Laboratory of Eco-Environment and Meteorology for Qinling Mountains and 
Loess Plateau, Xi' an 710014, Shaanxi, China; 2. Qinghai Climate Center, Xining 810001, Qinghai, China; 
3. Shanxi Province Atmospheric Observation and Technical Support Center, Xi an 710014, Shaanxi, China) 


Abstract: The prediction performance of four seasonal prediction model systems (BCC CSM 1.1, 
ECMWF System 5, CFSv 2, and TCC MRI-CGCM 3) was evaluated from a deterministic perspective. Focusing 
on the spatial distribution and temporal variation of precipitation climatology in early winter, BCC CSM 1.1 
exhibited the best prediction performance among the models. Additionally, TCC MRI- CGCM 3 performed well 
in capturing the interannual variability of precipitation, followed by BCC CSM 1.1. All models effectively 
simulated regional-uniform precipitation, with BCC CSM 1.1 exhibiting the highest time correlation coefficient 
and TCC MRI- CGCM 3 showing the highest pattern correlation coefficient. Although all models could 
reproduce the north-south reverse mode, BCC CSM 1.1 outperformed the other models in terms of reproducing 
the spatial pattern and interannual variation, given the model’ s strong prediction capability for the Eurasian (EU) 
teleconnection pattern and western Pacific subtropical high, as well as accurate simulations of the physical 
processes of El Niño- Southern Oscillation (ENSO) and Indian Ocean dipole (IOD) that impact early winter 
Tibetan Plateau precipitation. 
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